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In this communication with computer simulation we evaluate simple cubic random- 
site percolation thresholds for neighbourhoods including the nearest neighbours 
(NN), the next-nearest neighbours (2NN) and the next-next-nearest neighbours 
(3NN). Our estimations base on finite-size scaling analysis of the percolation prob- 
ability vs. site occupation probability plots. The Hoshen-Kopelman algorithm has 
been applied for cluster labelling. The calculated thresholds are 0.1372(1), 0.1420(1), 
0.0976(1), 0.1991(1), 0.1036(1), 0.2455(1) for (NN + 2NN), (NN + 3NN), (NN + 
2NN + 3NN), 2NN, (2NN + 3NN), 3NN neighbourhoods, respectively. In contrast 
to the results obtained for a square lattice the calculated percolation thresholds de- 
crease monotonically with the site coordination number z, at least for our inspected 
neighbourhoods. 
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1. Introduction 

Percolation [TJ5] is a mathematical description of a geometrical phase transition. This 
allows for purely theoretical studies of all phenomena occurring near the critical point 
with computer experiments solely or, sometimes, even analytically 3 . The mixture of 
occupied and empty sites of given lattice may exhibit some features of real physical 
systems. Among typical applications of the percolation theory one may find material 
science [J], immunology [5] or forest fires problems [B] and studies of liquids moving in 
porous media [7J, etc. [5J[H]. 

Generally speaking, the percolation theory deals with statistical properties of the 
clusters of occupied nodes (site percolation) or occupied edges (bond percolation) for a 
given graph, network or regular lattice. In the site percolation problem, the cluster is 
defined as a group of the occupied lattice vertexes which are direct or indirect neigh- 
bours. When each site is occupied with some probability p there is a critical probability 
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of sites occupation p c above which a cluster spanning through the whole system appears 
for the first time [HIS]- This special probability is called percolation threshold p c and 
it separates two phases (in the language of material science a conductor and an isola- 
tor). The value of percolation threshold p c depends on kind of percolation (site/bond), 
lattice/graph/network topology and assumed sites neighbourhoods. In the simplest case 
only the nearest neighbours constitute the neighbourhoods (von Neumann's neighbour- 
hood) or the nearest neighbours and next-nearest neighbours are considered (Moore's 
neighbourhood) . 

In the vicinity of the phase transition the quantity A describing the system follows a 
scaling relation 

A oc L a ■ f(x ■ L ), (1) 

where L describes the linear size of the system, /(•) is a scaling function and x is di- 
mensionless scaling parameter [TO] . For physical systems x usually plays the role of 
reduced temperature (T — Tc)/Tc, where Tq stands for critical temperature. The a 
and j3 parameters are universal exponents which — in the first approximation — do not de- 
pend on system details (kind of order/disorder phenomenon, lattice shape, site or bond 
percolation, etc.) but only on the system dimensionality [8,10 . However, for a and j3 
calculations the precise value of Tq is required. For the geometrical model of the phase 
transition percolation threshold p c plays the role of critical temperature. 

In this communication we evaluate with computer simulations |11) the random-site 
simple cubic percolation thresholds for neighbourhoods including the nearest-neighbours 
(NN), the second- nearest neighbours (2NN) and the third-nearest neighbours (3NN). Our 
estimations base on finite size scaling analysis [8lll0j of the percolation probability vs. 
site occupation probability plots. The Hoshen-Kopelman algorithm [T2] has been applied 
for cluster labelling. The calculated thresholds p c are 0.1372(1), 0.1420(1), 0.0976(1), 
0.1991(1), 0.1036(1), 0.2455(1) for (NN + 2NN), (NN + 3NN), (NN + 2NN + 3NN), 2NN, 
(2NN + 3NN), 3NN neighbourhoods, respectively. In contrast to the results obtained 
for a square lattice [13rH5] the calculated percolation thresholds decrease monotonically 
with the site coordination number z, at least for our inspected neighbourhoods. 



(NN) 




Fig. 1: Basic neighbourhoods for simple cubic lattice. The NN neighbourhood is often 
referred to as von Neumann's neighbourhood while combination of (NN+2NN) is called 
Moore's neighbourhood. We propose to name (NN+2NN-I-3NN) neighbourhood in simple 
cubic lattice the Rubik's neighbourhood. 
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2. Calculations 

We use Hoshen-Kopelman algorithm [TJ] for occupied sites labelling. In the Hoshen- 
Kopelman scheme each site has one label: all sites in a given cluster have the same labels 
and different clusters have assigned different labels. 

Examples of percolation probability P vs. sites occupation probability p for various 
neighbourhoods and for various lattice linear sizes L are presented in Fig. [5J We use 
finite-size scaling analysis to determine p c numerically. As for finite systems the phase 
transition is never sharp we observe it when for increasing lattice sizes L the P{p) curves 
become more and more steep and intersect at a common point corresponding also to the 
case of L — > oo 16J. The common cross-point predicts percolation threshold p c . Such 
strategy was successfully applied for many system description where phase transition 
may be observed including percolation [17] . Ising model [18] . majority- vote models [19] 
or opinion dynamics |20) . 

Basing on P(p) dependence for various L we look for an interval of the length Ap = 
2 ■ 10~ 4 where curves for L — 63 and 100 cross each other. The results have been 
averaged over TV = 10 5 lattice realisations. For example for 2NN case this interval is 
(0.1990,0.1992) and the centre of this bracket plays the role of percolation threshold. 
Basing on a uncertainty type B evaluation procedure |21j we estimate the uncertainty 
u(p c ) = Ap/V3 w lO" 4 . 

3. Results 

The P(p) dependencies for L = 22, 63 and 100 are presented in Fig. [2j The evaluated 
percolation thresholds p c for various neighbourhoods are collected in Tab. [U To check 
the accuracy of our estimations we evaluated the percolation threshold for compact neigh- 
bourhoods as well. The obtained values p c (NN) = 0.3116(1), p c (NN+2NN) = 0.1372(1) 
and p c (NN+2NN+3NN) = 0.0976(1) agree nicely with the results of extensive numerical 
simulations [35] and the earlier estimations |2"3"ll2~i] . 

Table 1: Simple cubic lattice random-site percolation thresholds p c for various neigh- 
bourhoods constructed with basic neighbourhoods NN, 2NN and 3NN. 



neighbourhood 


z 


Pc 


earlier estimations 


NN 


6 


0.3116(1) 


0.31160- •• [22] 


2NN 


12 


0.1991(1) 




3NN 


8 


0.2455(1) 




NN+2NN 


18 


0.1372(1) 


0.137 23 . 0.13735(5) [H] 


NN+3NN 


14 


0.1420(1) 




2NN+3NN 


20 


0.1036(1) 




NN+2NN+3NN 


26 


0.0976(1) 


0.097 [23], 0.0976445(10) [24] 




Fig. 2: Percolation probability P vs. site occupation probability p for various neighbour- 
hoods in three dimensions for three lattice sizes L. The results presented here have been 
averaged over N lattice realisations. 
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Fig. 3: Percolation thresholds p c vs. sites coordination number z for our (x) inspected 
three-dimensional neighbourhoods and other three-dimensional lattices [25-31] estima- 
tions (+). The straight lines are the least squares fits of p c {z) oc z -7 to the experimental 
data. 



These p c values decrease with sites coordination number z as shown in Fig. [3] The 
p c (z) dependencies may be fitted nicely by a straight line in a logarithmic plot, namely: 
p c (z) oc z -7 with 7 = 0.790(26). The percolation thresholds for some three-dimensional 
lattices have been used for fitting p c {z) dependence including ice [25], diamond [26II27] . 
hpc [28], bcc [23123130], fcc EUGl] and La 2 _^Sr :c Cu04 [31] lattices. These additional 
percolation threshold values are indicated as pluses (+) in Fig. [3j 



4. Conclusions 

In this communication for the first time we evaluate the random-site percolation 
thresholds for the simple cubic lattice with 2NN, 3NN, NN+3NN and 2NN+3NN neigh- 
bourhoods for which sites from the first, the second and the third coordination shells were 
included (see Tab. [1]). The obtained thresholds p c decrease monotonically with sites co- 
ordination number z according to the power law p c oc z~ 7 , with exponent 7 = 0.790(26). 

In contrast to the results obtained for a square lattice [TB] the calculated percola- 
tion thresholds decrease monotonically with the site coordination number z, at least for 
inspected neighbourhoods. 

The obtained results may be helpful in studies of the universal formulae [32] for 
percolation thresholds p c dependence on sites coordination number z. 

Finally, we propose to name (NN+2NN+3NN) neighbourhood in simple cubic lattice 
the Rubik's neighbourhood as it is identical with the famous Rubik's cube [33j — a very 
popular logical puzzle in early 80's. 
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